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ABSTRACT 

We compare the dark matter halos’ structural parameters derived for four Milky Way dwarf 
spheroidal galaxies to those of subhalos found in cosmological iV-body simulations. We confirm that 
estimates of the mass at a single fixed radius are fully consistent with the observations. However, when 
a second structural parameter such as the logarithmic slope of the dark halo density profile measured 
close to the half-light radius is included in the comparison, we find little to no overlap between the 
satellites and the subhalos. Typically the right mass subhalos have steeper profiles at these radii 
than measurements of the dSph suggest. Using energy arguments we explore if it is possible to solve 
this discrepancy by invoking baryonic effects. Assuming that feedback from supernovae can lead to 
a reshaping of the halos, we compute the required efficiency and find entirely plausible values for a 
significant fraction of the subhalos and even as low as 0.1%. This implies that care must be taken not 
to exaggerate the effect of supernovae feedback as this could make the halos too shallow. These results 
could be used to calibrate and possibly constrain feedback recipes in hydrodynamical simulations. 
Subject headings: galaxies: dwarf - galaxies: kinematics and dynamics 


1. INTRODUCTION 

The ACDM currently favored paradigm of structure 
formation has proven to be a successful model to explain 
the evolution and large scale structure of the Universe. 
It is on small scales, however, where the concordance cos¬ 
mology faces most of its challenges. One example is the 
well-known mismatch in the number of observed satel¬ 
lites around systems like the Milky Way and those found 
in cosmological iV-body simulations (Klypin et al. 1999; 
Moore et al. 1999). Although this can be solved by in¬ 
voking baryonic effects, a new problem has been reported 
regarding the abundance of massive satellites (the so- 
called “Too Big to Fail problem”, Boylan-Kolchin et al. 
2011). Pure Wbody simulations in ACDM also predict 
self-similar halos that have divergent central densities 
p ~ r -1 (Navarro et al. 1996b, 1997). However, mod¬ 
els of the internal kinematics of Milky Way dSph have so 
far yielded inconclusive results regarding the inner slopes 
of the halos hosting these galaxies (Battaglia et al. 2008; 
Walker et al. 2009; Breddels & Helmi 2013; Amorisco 
et al. 2014). 

The fundamental question that these studies have tried 
to address is whether a cusped dark matter density pro¬ 
file fits the observations better than a cored one. This 
is however, very difficult to answer with current data. 
For instance, Breddels et al. (2013) using non-parametric 
orbit-based dynamical models showed that, although 
small values of the logarithmic inner slope of the dark 
matter density profile are preferred in the case of the 
Sculptor dSph, the inner slope value is strongly degener¬ 
ate with the scale radius and the uncertainties arc large. 
A more feasible task is to measure the logarithmic slope 
at a finite radius, e.g. at r_ 3 , the radius where the loga- 
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rithmic slope of the light profile is —3. Breddels & Helmi 
(2013, hereafter BH13) have determined the logarithmic 
slope robustly (i.e. independently of the assumed profile, 
be it cored or cusped) at r_ 3 , for the dSph Fornax (Fnx) 
and Sculptor (Scl), and placed strong constraints on the 
values for Sextans (Sxt) and Carina (Car). In addition, 
BH13 confirmed using Schwarzschilcl modeling, that the 
mass at this radius can also be tightly constrained (Stri- 
gari et al. 2008; Walker et al. 2009; Wolf et al. 2010). 

This result shifts the discussion to a better posed ques¬ 
tion, namely: Are the logarithmic slope and mass at r _ 3 
inferred for these four dwarfs consistent with those found 
in simulations? This phrasing has the additional advan¬ 
tage that measuring these two quantities in simulations 
of dwarf galaxies is less prone to numerical artifacts. In 
this Letter we focus on answering this question using 
the Aquarius simulations (Springel et al. 2008). Since 
these simulations do not include baryons, discrepancies 
could possibly be attributed to missing physics. For in¬ 
stance, early adiabatic contraction could make the den¬ 
sity distribution of the central parts steeper (Blumenthal 
et al. 1986). On the other hand, feedback by supernovae 
(SNe) may modify the dark matter profile and lead to 
a shallower inner slope, especially for low mass galaxies 
(Navarro et al. 1996a; Mashchenko et al. 2006, 2008; Gov- 
ernato et al. 2012; Pontzen & Governato 2012; Teyssier 
et al. 2013; Di Cintio et al. 2014a, b). Finally, dynamical 
friction of large baryonic clumps has also been shown to 
modify the logarithmic slope in the same sense (El-Zant 
et al. 2001; Goerdt et al. 2006; Nipoti & Binney 2015). 

This Letter is structured as follows. In Section 2 we 
compare the observational results from BH13 to the sub¬ 
halos from the Aquarius simulations and highlight the 
discrepancies. In Section 3 we use analytic arguments 
to estimate the minimum amount of energy needed to 
transform the subhalos’ profiles to a distribution more 
consistent with the observational constraints following 
the arguments laid out by Penarrubia et al. (2012); 
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Amorisco et al. (2014); Maxwell et al. (2015), and derive 
the required SNe feedback efficiencies. Finally in Section 
4 we discuss the broader implications of our results. 

2. COMPARISON OF THE DARK HALOS STRUCTURAL 
PROPERTIES 

The Aquarius simulations are a suite of six (labeled 
Aq-A to Aq-F) high resolution cosmological dark mat¬ 
ter only runs of Milky Way-like halos. Here we use the 
highest resolution simulations available for all six halos 
(known as level-2), with a particle mass of ~ 1.4x 1O 4 M 0 
and a softening length of 66 pc (for more details see 
Springel et al. 2008). For a more direct comparison 
of the different systems, we have scaled all halos such 
that the Milky Way-like hosts have a common mass of 
Afhaio = 8.2 x 1O 11 M 0 (as in Vera-Ciro et al. 2013). 
A dark halo potentially hosting a galaxy such as Fnx 
(M v ; r ~ 1O 9 ' 5 M 0 ) is resolved with ~ 2 x 10 5 particles, 
and we expect our results to be reliable beyond a radius 
of ~ 200 pc, which for Fnx corresponds to ~ 0.2 r_ 3 . 

As mentioned in the Introduction, BH13 have shown 
that, besides M_ 3 (the mass within r_ 3 , Strigari et al. 
2008; Walker et al. 2009; Wolf et al. 2010), also 7-3 = 
dlog/ 9 dm/dlogr| r _ can be robustly determined for the 
dSph using current data (see also Wolf & Bullock 2012). 
7_3 does not depend strongly on the assumed paramet¬ 
ric profile that is used, whether cusped, cored or loga¬ 
rithmically cored 3 . These same two quantities can also 
be estimated for our simulated halos. For each Aquar¬ 
ius satellite we compute both M _ 3 and 7-3 using the 
values of r _ 3 measured for each individual dSph galaxy. 
M _3 follows trivially from the number of particles found 
within r_ 3 . We determine the slope 7 numerically using 
the algorithm described in (Lindner et al. 2015). The val¬ 
ues obtained are in good agreement with those derived 
assuming parametric forms for the density profile (such 
as truncated Einasto or NFW profiles), which are fitted to 
the data and from which 7 can then be computed. Com¬ 
paring the non-parametric and parametric estimates we 
gauge the uncertainty on 7 to be ~ 0 . 2 . 

In Fig. 2, we show the posteriors for the joint distri¬ 
bution of M _ 3 and 7-3 derived by BH13 separately for 
each dSph galaxy. The red contours correspond to the 1 
and 2 <j equivalent contours for the NFW (Navarro et al. 
1996b, 1997) profile 

Pdm(r) = P0 r/r s (l + r/r s )*’ (1) 

while green corresponds to the logarithmic cored profile 
(core 13 in the nomenclature of BH13) 

Pdm(^) PoTi i T~ M3 ’ (^) 

(1 + r/r () 3 

where r s and r' s are the scale radii and p d and p'o 
characteristic densities. We also plot in these pan¬ 
els the subhalos’ structural parameters determined non- 
parametrically and computed at the r _ 3 of the corre¬ 
sponding dSph. We show here only those subhalos that 
host a luminous galaxy according to the semi-analytic 

3 A distinction is made between strictly cored (lim r _»o dp 3m /dr 
= 0), and logarithmically cored (lim r _m d logPdm/dlog r = 0). 


model of Starkenburg et al. (2013). The black dots cor¬ 
respond to all such subhalos present in the six scaled 
Aquarius simulations, while the colored triangles are 
those that host a galaxy with an absolute magnitude 
within ±0.75 of the respective dSph. 

From the top left panel, we see that the subhalos over¬ 
lap little with Fnx, especially if only the correct luminos¬ 
ity range is considered. Note that if one only compares 
the M _ 3 values, then some subhalos with the right lu¬ 
minosity do actually match. Also if only the slope 7-3 
were to be taken into account, several subhalos would be 
compatible with the observational constraints. 

However the discrepancy becomes noteworthy for Fnx 
when comparing the observations and simulations in the 
two-dimensional space spanned by M _ 3 and 7 _ 3 . For Scl 
and Car there is some, although small, amount of overlap 
with the subhalos, but this is less so for Sxt. Both Car 
and Sxt have large uncertainties on the slope estimate, 
which makes the comparison less conclusive. 

Therefore it is clear that these four dwarfs all show 
a systematic deviation from the simulations: the slope 
of their dark matter profiles are shallower than those 
found in the Aquarius simulations, even if subhalos with 
a broader range of luminosities are considered. 


3. RESHAPING DARK MATTER HALOS 

As we have just seen dark matter only simulations in 
the ACDM cosmology, predict steeper slopes of the dark 
matter profiles than those determined for the dSphs. The 
most commonly proposed mechanism to transform the 
shape of the dark matter profile is feedback from SNe 
that affects the baryons and, through gravity, ultimately 
alters the dark matter mass distribution in these systems 
(Mashchenko et al. 2006, 2008; Governato et al. 2012; 
Pontzen & Governato 2012; Teyssier et al. 2013; Di Cintio 
et al. 2014b). 

To establish whether a given Aquarius subhalo can be 
reshaped into a halo that is compatible with a dSph, we 
first need to know how much energy is required to do the 
transformation. Following Peharrubia et al. (2012) we 
assume that each dwarf is in virial equilibrium, and we 
compute the energy needed to transform a subhalo i to 
match the dSph model j (given by Eq. (1) and Eq. (2)): 
SE = (WdSph — ffsub) /2. Here W su b is the potential 
energy of the subhalo and MyjSph is the potential energy 
of the best fitting model for our dSph: 

pT'ma, x,i 

WdSph = - 47 tG / drp dm ^(r)Mj(r)r, (3) 
Jo 

with M(r) = Md m (r), i.e. we ignore the baryonic com¬ 
ponent which is sub-dominant for dSph. We choose two 
different values of r max , namely r$o and rgs, the radii en¬ 
closing 50% or 95% of the bound mass of each subhalo. 
The potential energy at rg$ is more negative than at rso, 
and thus the energy required to transform the system in¬ 
side r 50 is smaller than that required for rg 3 . (Note that 
these values are larger than those used by Maxwell et al. 
(2015), who took r max = 3r c , with r c the radius where 
the dark halo profile’s logarithmic slope is —0.5). 

We calculate W su b directly using the particles from the 
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Fig. 1.— Joint probability distribution function for the logarithmic slope 7-3 and mass M— 3 of the dark matter halo at r _3 for our 
sample of dSph. The contours represent the lcr and 2<j-equivalent confidence regions for an NFW profile (red, Eq. (1)) and the core 13 model 
(green, Eq. (2)). The black points correspond to the subhalos of the six scaled Aquarius Milky Way-like halos, while the colored triangles 
are the subset with luminosity within ±0.75 of the value of the respective dSph as predicted by the semi-analytic model of Starkenburg 
et al. (2013). The error bar on the bottom right indicates the uncertainty in the slope estimated for the simulations. 


dark matter simulations: 

r'T max,i 

Wsub = — 47 tG / drp dm! j(r)Mj(r)r 

Jo 

pT max,! 

= —G I d MMi(r)/r 
Jo 

N 

= -^y> 2 k/r k , (4) 

k =1 

where m is the dark matter particle mass, rk is the sorted 
radial distance of particle k from the center of the sub¬ 
halo such that tt is the innermost particle and r jy that 
farthest away but still within r max . 

Fig. 2 shows for each dSph the ratio WyiSph/Wsub as a 
function of absolute magnitude My as given by the semi- 
analytic model. The blue bands indicate a magnitude 
range of ±0.75 around the corresponding dSph. The solid 
circles are for the NFW (red) and corel3 (green) best fit 
models computed for rgs while the triangles are for r$ o- 
When the ratio ll'dSph/IFsub < 1 (all symbols below 
the black horizontal line) objects need energy to trans¬ 
form their halos. In these cases the binding energy of the 


resulting object is larger (less negative) than that of the 
original untransformed subhalo. When WdSph/Wsub > 1 
(all symbols above the black horizontal line) the binding 
energy of the resulting object is smaller and such objects 
would have to lose energy to be transformed into a profile 
that matches the observations. Regardless of the mech¬ 
anism, this transformation is energetically possible and 
we do not discuss it further. 

As in Penarrubia et al. (2012), we estimate the total 
available energy from SNe as 

M* 

A’sn, total = 7 -c£(m* > 8 M q )I?sn, (5) 

{m*} 

where £(m* > 8M 0 ) = 0.0037 and (m*) = O.4M 0 . If the 
energy needed to transform a given halo is larger than 
what can be delivered from SNe (i.e. the SN efficiency 
is larger than 100%), we plot it in Fig. 2 with an open 
symbol. Note that for Fornax only a few open symbols 
exist. This means that almost all of the candidate subha¬ 
los could be reshaped into Fornax given its stellar content 
and the energy available from SNe. 

For each dSph, we plot in Fig. 3 the ratio of the en¬ 
ergy required to transform the system to the SN energy 
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Fig. 2. — Ratio of the potential energy of the dSph, WdSph? t° 
that of the Aquarius subhalos Wg U b> as a function of their pre¬ 
dicted absolute magnitude My, computed for 7*95 (circles) and 7*50 
(triangles) for the NFW (red) and core 13 (green) models. The blue 
band indicates an absolute magnitude range within d=0.75 of the 
of the corresponding dSph. The open symbols are subhalos which 
cannot be transformed with the energy available in SNe. 

output as estimated above, i.e. esn = 52£/.Esn, total is 
the SN efficiency needed to do the transformation. We 
thus only plot the objects which lie below the horizontal 
lines in Fig. 2. This analysis reveals that for each dSph, 
sublralos with the right luminosities exist which can be 



Fig. 3.— Similar to Fig 2, now showing the SN efficiency needed 
to transform the Aquarius halos. 

reshaped to follow a dark matter profile consistent with 
the observations. The required efficiencies range from 
~ 0 . 1 % up to slightly more than 100 %, depending some¬ 
what on the extent of the region of the subhalo that 
would be transformed (rso or ^ 5 ). Furthermore there is 
always at least one subhalo that requires a SN feedback 
efficiency £sn < 0 . 001 . 

4. DISCUSSION AND CONCLUSIONS 
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We have shown that in the two dimensional parameter 
space defined by the mass and the logarithmic slope of 
the dark matter density profile at r_ 3 , the Milky Way’s 
dSph do not match the predictions from dark matter only 
simulations. While in Vera-Ciro et al. (2013) we argued 
that by rescaling the Milky Way mass to a (low) value of 
~ 8x 10 11 Mq we could circumvent the so called “Too Big 
Too Fail Problem”, we now see that halos with the right 
mass, do not have the right logarithmic slope at r_ 3 . 
However, baryons can have a considerable effect on the 
shape of the dark matter density profile through feedback 
(e.g. Governato et al. 2012). We have shown here that it 
is energetically possible at redshift of z = 0 to reshape 
the simulated halos and make them compatible with the 
observations. Note that this process would require even 
less energy at higher redshift, as discussed by Amorisco 
et al. (2014); Madau et al. (2014). 

Hydrodynamical (cosmological) simulations of (iso¬ 
lated) dwarf galaxies are now available in the literature 
(e.g. Teyssier et al. 2013; Di Cintio et al. 2014a, b; Sawala 
et al. 2015), and are more representative of galaxy for¬ 
mation scenarios in ACDM, since they self consistently 
include baryonic physics. Di Cintio et al. (2014b), for 
instance, predicts the slope in a given radial range to be 
predominantly dependent of the ratio M*./Mhaio- In their 
simulations an object like Fnx (whose M* = 4x 10 7 Mq, 
de Boer et al. 2012), is embedded in dark matter halo of 
virial mass Mhaio ~ 10 10 — 10 lo - 5 M Q . Such objects are 
reported to have dark halo profile’s logarithmic slopes 
a £ [—1, —0.5] for the radial range 0.8 — 1.8 kpc. This is 
inconsistent at a 2 er level with our measurement of the 
slope for Fnx at r _3 ~ 1 kpc as can be seen from Fig. 2. 
Such simulations therefore, predict slopes that are too 
shallow compared to the observationally determined val¬ 
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ues (see also Madau et al. 2014; Chan et al. 2015). This 
suggests that the feedback scheme used might be too 
strong for objects of this mass scale. 

It is more difficult to make a detailed and conclu¬ 
sive comparison the other dSph because the simulations’ 
suites available in the literature do not contain such faint 
objects or the ranges in which the slopes are measured 
are typically larger than the observed r_ 3 . 

The steep slopes (corresponding to high concentration) 
of low mass dark halos could be another manifestation of 
the “Too Big to Fail” problem (however see Papastergis 
et al. 2015), although the two issues are not necessarily 
strictly related, since other solutions to this problem 
exist (Wang et al. 2012; Vera-Ciro et al. 2013; Brooks 
& Zolotov 2014; Arraki et al. 2014). For example, 
Brooks & Zolotov (2014) argue that efficient feedback 
solves it, but the preliminary comparisons made above 
show that a careful treatment is required. We hope 
that our measurements of the logarithmic slopes and 
masses of dSph can be used to better constrain poorly 
understood processes such as the interplay between 
adiabatic contraction and supernova feedback in these 
small mass systems. 
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